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Abstract 

Molecular clouds in the interstellar medium suffer gravitational instabilities that 
lead to the formation of one or multiple stars. A recently formed star inside a cold 
cloud communicates its gravitational force to the surrounding environment and soon an 
accretion flow falling into the star develops. After their formation, all stars soon eject 
a wind of gas that interacts with the external accretion flow. This interaction produces 
a shock wave that evolves with time. The work presented in this article formulates a 
simple prescription for the evolution of this interaction. With the aid of this model we 
construct a few radio continuum maps of the source. 
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11 Introduction 



When a low mass, isolated star is formed inside a gas cloud because of a gravitationa l 



collapse, the process of accreting gas from the cloud onto the star begins ()Shu et al 



19871) 



At a later time a wind of gas emanating from the star collides with the accretion flow. 
Depending on the ram pressure that the wind generates, it will be able to expand forever 
(when the ram pressure of the wind is much greater than the ram pressure of the medium) , 
reach a steady state (when accretion and wind ram pressures are comparable) or collapse 
in the surface of the star (when the accretion ram pressure is greater than that of the wind) 



(Wilkin fc 



Stah 



Bondil 



er 



199^ . 



2003). 

1952t ) investigated the form of a spherical accretion flow induced by the grav- 
itational pull of a central object. Because of the spherical symmetry of the accretion 
quantities, the shocks produced by th e sphe rical wind and the accretion flow are both 



spherical and evolve with time. 



UlrichI (|1976r ) made modifications to the spherical accre- 



tion model by considering a cloud with a small angular momentum. The cloud is rotating 
about an axis on which the central star lies. The result is that the accretion quantities in 
this model have cylindrical symmetry, i.e. the y do not depend explicitly on the azimuthal 



angle. The interaction between this ' 



UlrichI flow" and the spherical wind from the star 



produces two shocks that evolve with time. In this article we develop a simple model for 
this particular hydrodynamical interaction. We do this by equating the ram pressure of 
the wind with the ram pressure of the accretion flow. This is the easiest way to model the 
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interaction. 



Wilkin &: Stahlen (|2003) have already constructed a numerical solution to the problem 
by solving the complete hydro dynamical equations. These authors have included on their 
code not only pressure forces that give rise to the shocks. On their analysis they also 
included centrifugal forces that appear because of the curvature of the shock, and they 
took into account the gravitational field produced by the central star. The problem with 
such a general solution is that it is not very easy to predict observational quantities. In 
the present article we show how, with a simplified method, it is possible to produce radio 
continuum maps of the interaction between the wind and the accretion flow. 



In section §2 we describe briefly 



Ulrich's accretion model. Sections S3 and §4 describe 



the steady state and the time evolution of the interaction. Finally, in section §5 we describe 
the validity of the predicted thin layer approximation and we calculate radio continuum 
maps from our model. 



\2 Ulrich's accretion model 



Let us briefly discuss some of the most important properties of a small angular momentum 
accretion flow that falls to a central c ompa ct object, such as a star. This type of accretion 



flow was first investigated by 



UlrichI Hmd ) 



Consider a star, or any other compact object of mass M. The star is embedded on a gas 
cloud, which for simplicity will be consider to be of infinite extension. Far away from the 
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star, the density poo, the pressure poo and the velocity of sound Coo have uniform constant 
values. We assume that the cloud is rotating about the z axis as a solid body in such a way 
that its specific angular momentum F has the value T ^ far away from this axis. We will 
consider the case in which the accretion process is so low, that the mass M of the central 
object can be taken as constant. Finally, we assume that the flow has reached a steady 
state. The accretion flow obeys a polytropic relation 



where p and p are the pressure and density of the gas respectively, and k is the polytropic 
index. 



and under the assumption that the mass of the central object is much greater than the 
gas mass contained on a sphere of radius r (r the distance from the star to a certain fluid 
particle), the self gravity of the gas is negligible with respect to the gravity of the central 
object. In what follows we will consider Fqo to be small so that the accretion with rotation 
can be considered to be a perturbation from the non-rotational accretion Bondi flow. 

The problem is characterised by the following parameters: the gravitational constant G, 
the radius of the central object ro, the constants Fqo, M, poo, Poo, Coo, and the polytropic 
index k. With these quantities it is possible to build three dimensionless parameters 




(2.1) 



Due to the fact that the accretion rate M in iBondil (jl^^^) accretion flow is constant, 
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(without counting k) 



(2.2) 



GM' V GM 

and Poo/PooC^ ~ 1- The length Vc = GM/c^ is (apart from a function that depends of k 
which is of order one) the radial distan ce whe r e the fluid particle reaches the velocity of 



sound in the accretion without rotation (|Bondi , 



19521 ). In other words, the fact that 5 <C 1 



ensures that the radius of the central object is small compared with this distance. 

The dimensionless parameter e is ratio of the centrifugal force to the gravitational force 
of a fluid particle evaluated at position re- Thus, if we impose the condition e ^ 1, the 
rotational effects are small perturbations to the Bondi accretion flow. 

Due to the fact that the angular momentum of the accretion flow must be conserved, 
eventually the centrifugal force (~ F^/r^) will be equal to the gravitational force (GM/r^). 
The position = ere where this occurs is precisely the distance where the flow deviates 
considerably from the Bondi spher i cal ac cretion and where a disc of radius on the 



equatorial plane forms (jCanto et al 



ma) 



Low mass, recently formed stars are such that tq ~ 2Rq, c^o ~ lkms~^, M ~ IMq, 
with an azimuthal velocity far away from the star of the order of 10~^^s~^, and a 
lifetime r value such that lO^yrs ^ t < lO^yrs. With these values, the angular momentum 
far away from the rotation axis takes the value Loo = (cr)^^? ~ 9.61 x 10^^ cm^ s~^. Thus, 
for low mass, recently formed stars, it follows that 5 ~ 1.75 x 10^^ ^ 1, e ~ 0.2 <C 1, 
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rc ~ 536 AU and ra ~ 100 AU. 

At time r = 0, when the gravitational effects produced by the star begin to be noticed 
by the gas in the cloud, the gas starts to be accreted. In what follows we will consider 
that the time r is sufficiently large so that Jq" dr/c{r) < r, in which c(r) represents the 
local sound speed. For instance, if the accretion gas is isothermal, the speed of sound is 
constant and the previous inequality means that Vc < ct. For the general case in which 
the velocity of sound is not constant, this argument implies that the points r contained 
inside the sphere determined by /([dr/c(r) > r are not perturbed by the gravitational 
effects of the star. On the other hand, the gas particle points r inside the shell that 
satisfies J^" dr/c{r) < dr/c{r) < r are such that they are being accreted to the centre 
subsonically and its rotational effects are negligible. The flow inside the region < r < Tc 
is being accreted with a supersonic velocity and its rotational effects are small. When 
r < Td, the rotational effects are very important on the dynamics of the accretion flow and 
the velocity of the gas flowing towards the centre is supersonic. 

Because pressure gradients and variations on the internal energy along the streamlines 
of the accretion flow inside the sphere of radius rc contribute little to the energy and 
momentum equations, the streamlines can be taken as ballistic. If we also consider that 
the mass of the disc is much smaller than the mass of the central object, then the 
trajectory of a fluid element of the accretion flow is described by a central Newtonian 
potential. 

The total mechanical and internal energy of a fluid element far away from the central 
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object is non-zero. This is because the azimuthal component of the velocity and the 
temperature differ from zero at that point. Nevertheless, the difference between the total 
energy and the kinetic plus potential energy of a fluid element near the disc is small. This 
is valid for small values of the rotational kinetic energy (i.e. when e <^ 1) and when heating 
due to radiation is insignificant. In other words, the streamlines near the disc have to be 
parabolic trajectories such that 

E=-v^ «i 0, (2.3) 

2 r ^ ' 

where E is the energy per unit mass of a fluid element near the disc and v its velocity. Once 

the particles reach the plane ^ = 7r/2 {9 the polar angle) where the disc lies, the particles 

merge into the disc. Because of this, the only interesting trajectories for the accretion flow 

are the upstream trajectories to the disc. 

For simplicity, in what follows we will work with dimensionless quantities, unless stated 

otherwise. In other words, let us make the changes 
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where the natural velocity f^* and density po are defined by 



M 



The analytical form of the s tream lines, t he velocity field and density for the accretion flow 



with rotation are given by (jUlrich 



1973) 



sin^ 6*0 



1 — cos Q I cos 6*0 ' 



(2.6) 



1 + TTTTT ' (2-7) 



r J \ cos 6*0 

..=(^y'' ^°^'°."r' fi+^V'', (2.8) 

r J sm t) \ cos t^o / 

1 y/V^ _ «^y/^ sin^ (2 9) 

r) V cos 00 y sine' ^^'^^ 



p = r-3/2 I 1 + ^22:^ ) {l + 2r-^P2(coseo)} \ (2.10) 
V cos6'oy 

where Vi {i = r, 9, ip), and p are the velocity and density fields of the accretion flow respec- 
tively. The azimuthal angle is represented by ip. The initial polar angle Oq that a fluid 
element has once it begins its way down towards the disc (see Figure^, M is the accretion 
rate and P2 (^) is the second order Legendre polynomial given by 



n \ -1/2 

1 + ^ {l + 2r-'P2icos9o)}-' 

cos uo J 



P2{o = l{^e-i)- (2.11) 



*iik is the velocity (keplerian velocity) of a particle that rotates about an axis on a circular trajectory of 
radius rd under the influence of a central Newtonian field. 
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A plot of the streamlines (Eq. (|2.Hj) ) projected on a plane 99 = const is shown in Figure|21 
The angle labels a particular streamline. For each there is only one streamline, because 
they cannot intersect. It is clear from the velocity field equations that when ^0 ~^ 0, the 
polar and azimuthal components of the velocity become zero as well. Thus, the streamlines 
become parallel near the axis of rotation. 

The velocity and density field equations are functions that depend only on the particle's 
position and not of the initial polar angle ^o- This is easily seen by rewriting Eq. ()2.3|) as 



cos' 



;^ 6*0 + cos 9o{r - 1) - rcosO = 



(2.12) 



which has a solution (see appendix §6 1 




for r = 1 




for r > 1 



cos^o = < 



(2.13) 




for r < 1 and (l£|i^)2 _ [l^f > 




for r < 1 and ( 



r cos 6 \ 
2 / 



)'-(¥)'<o- 
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Thus, Eqs. (|2.7j) - (|2.1()|) and Eq. (|2.i:-{j) provide the velocity and density fields as functions 
of the polar angle 9 and the radial coordinate r only. 
From Eq. (|2.12() it is easy to see that 



cos 9q 



G(l -r)^/r 



cos ( 



1 ■ 

2 \r + 2 



as 



(2.14) 
(2.15) 



where 0(x) is the Heaviside step function with a value of 1 for x > ^-iid for x < 0. 
With the substitution of Eq. ^(TT^ and Eq. (fTT^ into Eq. (|m7|) it is possible to give 
values for the density in the equatorial plane and the polar axis respectively 



p{e = 7:/2) = { 



\r (1 - r) ^ , for r < 1, 

(2r - 1)"^/^ (r - 1)"^ , for r > 1, 
p(^ = 0) = (2r)-i/2 (r + 2)~^ 



(2.16) 



(2.17) 



The fact that the density tends to infinity as r tends to zero for any 9 is due to the 
accumulation of the accreted gas around the central object. Nevertheless, the border of 
the disc (i.e. r = 1 and 9 = vr/2) has also an infinite density. This is because we are 
considering an infinitely thin disc and so, border effects on it are expected to appear. 

Figure shows the variation of the density as a function of position for different polar 
angles. This plot, together with the density contours (as shown by Figure 0J) can be used 
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to make a detailed analysis of the density distribution in the accretion flow. 

§3 Steady interaction 

From now on, let us suppose that the central object is a recently formed star, and that a 
supersonic, spherical wind with constant radial velocity is being produced by it. In this 
section we discuss the interaction between the Ulrich accretion flow described in Section 



2 with the wind of the star. 

The interaction between the supersonic accretion flow and the supersonic stellar wind 
gives rise to an initial discontinuity. The end result of this interaction is the production of 
two shock waves and a tangential discontinuity between both. In the frame of reference of 
the tangential discontinuity, the shocks separate from each other. Under the assumption 
that both shocks occupy the same position in space, so that we can call both shocks 
"the shock" for simplicity, we analyse the geometrical form that the shock has under the 
assumption that the only type of forces that keep the shock stable are pressure forces. 
Apart from pressure forces, there will be centrifugal forces (produced by the post-shock 
material) acting on the walls of the shocks. These forces are originated because the shock 



is curved. In section §5 we analyse under which conditions, the assumption that both 
shocks occupy the same position in space (thin layer approximation) is valid by studying 
the cooling lengths of the flow. It turns out that for the steady case, this approximation is 
correct. 
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Under these assumptions, the locus of the points of the shock are given by balancing 
the ram pressure produced by both flows, the accretion flow and the stellar wind. In other 
words, the shock is located at points for which 

Pw^^wn = pvl, (3.1) 

where the subindex w labels quantities related to the stellar wind and n means the normal 
component of the velocity. The wind's density is determined by the mass loss rate of 
the star 

Mw = 47rr2vwPw (3.2) 
In order to get a dimensionless form of Eq. (|3.1|) . we make the changes 

V-w PwQ 

where = M^/Airr'^v^. Combining these changes with the ones presented in Eq. (|2.4p 
we obtain 

Pw^wn^ = P'^l (3-4) 
where A is a dimensionless parameter given by 



(3.5) 
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Eg . (13.41) determines the position r{9,(p) of the shock as a function of the polar angle 6 
and the azimuthal angle if. A normal vector to this surface is given by 

dr dr 

since dr/d9 and dr/dip are tangent vectors to any surface parametrised by the polar 
and azimuthal angles. Because the accretion flow and the stellar wind have azimuthal 
symmetry, the shock wave position should not depend on the angle (p. Since the length 
element dr = drer + rdOeg + r sin6'd99e^, where Cj (i = r, 9, ip) are unitary vectors in the 
i-th direction, then Eq. ()3.6|) can be rewritten in the following form 



{r2 + (ar/a0)2} 
where n = n/n is a unit vector in the direction of n 



The accretion flow and wind velocities are radial on the rotation axis. This means that 
the locus of the shock satisfies the following boundary condition 



dr 
89 



= 0. (3.8) 

9=0 



Using Eq. 1)3. 7() and Eq. (|3.8() . it follows from Eq. i)3.4|) that over the rotation axis of the 
cloud, the value of the distance r from the star to the shock is given by 



r 



(3.9) 
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This equation implies that for A > 1/2 there is no value for the initial condition of the 
position of the shock. In other words, there is no steady solution for the case when A > 1/2. 

The equation satisfied by the spatial points on the shock r{6) is given by substitution 
of Eq. (jSH) into Eq. ((231), that is 

d9 - • ^^-^"^ 

Figure shows integrals of Eq. 1)3. 1U() for different values of the parameter A. These 
integrals were obtained with a 4th rank Runge-Kutta method. From the figure it follows 
that when A < 1/2 the shocks can be classified in two cases. Those that have positive 
derivatives with respect to the polar angle (A = 0.2, for example) and those for which this 
derivative is negative (A = 0.48, for example). Let us show that there is no configuration 
for which this derivative is zero for all polar angles. For this, it is sufficient to assume that 
the configuration is a circumference. It then follows that r{9 = 0) = r{9 = Tr/2). From 
this and Eq. 1)3. 9|) it follows that the obtained value for A does not correspond to a circle. 

Until now, we have considered that both shocks, the accretion shock and the wind shock, 
occupy the same position in space. However, as shown by the cartoon of Figure IHl there 
is an intermediate region of shocked gas between them. Let us calculate the direction that 
the post-shocks should have immediately after crossing their corresponding shocks. Under 
the assumption that both shocks are highly radiative, the post-shock normal velocities are 
negligible. Nevertheless, since the tangential components of the velocity are continuous 
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through the shock, the post-shock direction is determined by the tangential components 
of the pre-shock velocity. Because the shock is not a function of the azimuthal angle, a 
unitary tangent vector to it's surface is given by 



Due to the symmetry of the problem, it is useful to analyse only the region < 9 < 7r/2. 
Using the tangent vector of Eq. H3.11() and because the polar component of the velocity is 
positive for the accretion flow, we say that the post-shock flow (either accretion or stellar 
wind) "goes up" if • r < (i.e., the flow leaves the equatorial plane) and the flow "goes 
down" if i> • T > (i.e., the flow approaches the equatorial plane). Figure Q shows values 
of the pre-shock tangential velocities for both flows as a function of the polar angle for 
different values of the parameter A. From this figure it follows that in the case of the stellar 
wind, for A < 0.3 the post-shock flow goes down and for A > 0.3 this post-shock flow goes 
up. When A ~ 0.3 the post-shock flow goes up in some regions and goes down in some 
others. For the accretion flow, the situation is different. It always goes down, regardless 
of the configuration chosen. It is important to note that there is an accumulation point of 
material, corresponding to the polar axis (^ = 0). In this point, both post-shock flows have 
no motion. In this accumulation point, it is expected that the post-shock flow collimates 
to regions away from the star when the centrifugal forces are taken into account. 



{dr/d9) Br + reg 



(3.11) 



T = 
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§4 Time evolution 

Let us now analyse the evolution in time of the interaction between the stellar wind and 
the accretion flow. As discussed in section [|3| two shocks are formed. For simplicity we 
will consider that both shocks occupy the same position in space and we will refer to this 
region as the shock, unless stated otherwise. If we assumme that only ram pressure forces 
act on the surface of the shock, then from the system of reference of the shock Eq. p. If) 
is also valid. In order to transform it to the frame of reference of the star, we make the 
transformation 

> {Vsn - Vn) , (4.1) 

where Uwn, Vsn, represent the velocities of the stellar wind, the shock wave and the 
accretion flow respectively, all normal to the shock wave. With these changes, Eq. (|3.1I) 
gives 

Pw {Vwn - Vsnf = P (Vsn " ^^n)^ • (4.2) 

In order to write Eq. 1)4.2(1 in a dimensionless form, we use the changes described in 
Eq. ((131) and Eq. (|^ together with 
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Performing these transformations, Eq. (|4.2j) can be written as 

where A is defined by Eq. and rj is another dimensionless parameter given by 



(4.4) 



When Vsn — *■ 0, that is when the evolutionary shock becomes steady, Eq. (|4.4)1 converges 
to Eq. ((231). 

Since we are only interested in the geometrical shape of the shock and not on the 
particular trajectory that a small fragment of the shock follows, we analyse only the form 
that the shock will have as a function of time for a fixed polar angle. Because of this, the 
relation between the shock's radial velocity and its normal velocity can be obtained from 
Eq. ipTTfl 



TVs 



\r'^ + {dr/def\ 
With this, Eq. 1)4. 4() can be rewritten as 



1/2- 



(4.6) 



dr _ (Apw)^''^ r + P^^"^ {rv^ - vedr/dO] 



(4.7) 



r 

where the normal values of the stellar wind and the accretion flow given by Eq. 1)3. 7|) have 
been used. 
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Near the star, for r <^ 1, the velocity field and density for the accretion flow are 
spherically symmetric and coincide with the spherical symmetry of these quantities for the 
stellar wind. In other words, at time t = the shock is a spherical surface of radius r^, the 
radius of the star 



r{0,t = O) = — < 1. 



(4.8) 



Since Eq. 1)4.4(1 has two dimensionless parameters, we will co nsider in what follows tha.t 



1983) 



the star is a recently formed star of low mass. For these stars ([Black fc Matthews 
fw ~ 300 km s""*^, M ~ 1M0 y ~ 100 au. With these values, the parameter r] can be 
determined uniquely by the values of A through the following expression 



(4.9) 



The solution of Eq. ()4.7 1 with the initial condit i on giv en by Eq. (|4.8|) was obtained by 



numerical integration using 



MaCormack &: Paullavl (|1972l ) corrector method. The integrals 



for this equation are plotted in Figure |HJ The solutions show that the steady state is only 
reached for A < 1/2, and so there is a convergence of these solutions to the ones obtained in 
the previous section. When A > 1/2 the shock grows indefinitely, apart from the equatorial 
plane, where the shock reaches the border of the disc and stays there for sufficiently large 
times (see the case A = 0.7 of Figure |H1 for example). This effect is due to the infinite value 
that the accretion density reaches in the border of the disc. 
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Using Eqs. (|2.Hj) - (|2.1()|) for r S> 1 and Eq. (|4.7j) it is possible to obtain the relation that 
the shock surface must have for sufficiently large times 



which has a solution 



dr 21/4A1/2 



dt 



.1/4 



as r ^ oo, 



(4.10) 



5AV2 

27/4 



4/5 



(4.11) 



Figure IHl shows for which times the approximation made with Eq. H4.11() is valid on the 
polar axis of the shock. 



§5 Astrophysical consequences 

We have assumed that the shock produced by the stellar wind, as well as the one produced 
by the accretion flow, occupy the same position in space. In this section we analyse under 
which circumstances this approximation is correct. As usual, we define the cooling length 
Ic as the distance travelled by a particle fro m the shock front un til it reaches a point in 



wh ich the tempera t ure h as a value ~ 10^ K ( Hartigan et al 



Hartigan et al 



19871 ) 



(1983) estimated cooling lengths for different shock velocities Vsn in the 



interval 20 km s ^ < < 400 km s ^, and particle number density n in the in terval 



100 cm -3 < n < 1000 cm . From their results it follows that (ICanto et al 



1983) 
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where 



and 



MJ ) ~ V lOOkms-i 



n 



cm 



-1 



(5.1) 



A= < 



6.6 Vsn < TOkms ^ 
43.33 > 70kms~^ 



(5.2) 



(5.3) 



-5.5 Vsn < 70kms ^ 
5.7 fsn ^ 70kms~"'^ 

Prom now on, we will assume that Eq. ()5.1|) is valid for any density and velocity ranges. 

The equations that describe the accretion flow (Eqs. (|2.6|) - H2.1U() ). the stellar wind flow 
(Eq. ()3.4|) ). the shock velocity (Eq. 1)4.4(1 ') and the time can be written in dimensional form 
as follows 



r = r^r, p = pop, Vi = V],Vi {i = r, 6, ip), = V]^Vs, 

v^ = v„v^, Pw = /OwOPw, t = {rd/v]^)t, (5.4) 

where tilded quantities are referred to the dimensionless variables considered in the previous 
sections. We now give values of the hydrodynamical quantities for the accretion flow and 
the stellar wind for typical values of recently formed stars of low mass 
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'0 

no \ a I M 

^ ^ 3.6 X 10^ ' 



cm"3; ■ \ 10-6MQyr-i 



rd \ -3/2 f M\~^ 



(5.6) 



( ) 
VIOOAU / 



Mo 



cm 



3 X loM TTT^Moyr- 



3/ I 10-7 



X 



(5.7) 



-2 / . N -1 



■) ( 



.100 AU / V 300km s-i- 
where no y rawo are the particle number densities of the accretion flow and the stellar wind 

respectively with the assumption that the average mass per particle is m ~ 1.3mi^. As 

shown by Eq. ()5.5() . typical values of the velocity fk are sufficiently small to be used in the 

calculation of the cooling length of the accretion flow using Eq. 1)5. 1|) . Thus, from now on, 

we will only analyse the cooling lengths for the stellar wind. 

In order to simplify things further we calculate the cooling lengths for the stellar wind 
leaving fixed the quantities M, M, r^, fw, with typical values as the ones showed in 
Eqs. H5.5|) ~ H5.7|) and we will only vary the quantity M^. Table shows the maximum 
variation of the ratio of the cooling length to the length between the shocks for different 
polar angles in the steady case. These results show that, for the steady case the assumption 
of a thin post-shock layer is correct for the stellar wind. 

Let us now analyse the cooling lengths of the stellar wind for the case in which the 
shock evolves with time. In the frame of reference of the shock wave, the formulae used 
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A 


Hc/r)^.. X 10=^ 


0.04 


0.54 


0.08 


1.08 


0.12 


1.65 


0.16 


2.25 


0.20 


2.88 


0.24 


3.56 


0.28 


4.27 


0.32 


5.20 


0.36 


6.38 


0.40 


7.96 


0.44 


10.35 


0.48 


15.21 



Table 5.1: Values of the cooling length in the steady case for the shock produced by 
the stellar wind as a function of the parameter A. The quantity {lc/r)niax represents the 
maximum value of the ratio of the cooling lengths Ic to the position of the shock r, for 
different polar angles. 
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for the steady case are valid. With the aid of Eq. (|4.1j) . in the frame of reference of the 
central star, Eq. ()5.1() takes the following form 

VAUy V lOOkms-i J \ cm-3 
In this case, it follows that the ratio of the cooling length to the radial distance to the 
star, when A > 1/2, is large. To see this, let us make an asymptotic expansion of Eq. ()5.8|) 
for r 1. With this, and using Eq. (|4.11jl . it is found that 

as i — > oo. Figure ITUl shows at which times the thin layer approximation ends its validity 
for certain values of the parameter A > 1/2. In the case where the shock evolves with time, 
for A < 1/2, the cooling lengths and the distance to the shock are always small until the 
steady configuration is obtained, converging to the values shown on Table 

Let us now analyse the shape that a radio continuum observation will show of the 
interaction described so far. For this, we assume that the shock emission is due to free-free 
processes and we use the brightness temperature Tb given by 




Tb = T {1 - e-^^) , (5.10) 

where T is the temperature produced by the shocks. The optical depth Ti, in the continuum, 
at a radio frequency of a shock wave with velocity Vs and a pre-shock density n is given 
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by (|Ciirie] et a1.l.ll989h 



T„ = 9.83 X 10" 



n 



cm 



ViO'^Ky 



300kms-i 

0.55 



5 GHz 



-2.1 



For our model, it follows that -^1 and so, Eq. H5.1U() can be written as 



(5.11) 



Tb = Ttu, as Ty 0. 



(5.12) 



Figures I11I12I show different projections of the brightness temperature isocontours on 
the plane of the sky for typical values of = 10^'^MQyr^"'^ (A = 10). 

Let us now calculate the emission flux produced by the stellar wind shock for the value 



used in Figures 1111121 in which Af„ = lO^'^MMyr "^ (A = 10 



flux received on the plane of the sky is given by (|Rodriguez 



For this case, the emission 



199(1 ) 



Sy = — — / dVtTyo, (5.13) 

where the integral is extended over the solid angle Q occupied by the shock surface on the 
plane of the sky and c is the speed of light. The optical depth t^q is observed through 
the plane of the sky. The quantity Ty^dD. = Tyda/D"^, where da is the area element and 
Ty is the optical depth, both normal to the shock surface. With this and because D is the 
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distance from the observer to the shock surface, Eq. can be rewritten as 



Using Eq. (|3.7j) in spherical coordinates, Eq. (|5.14|) takes the form 



S„ = f|^/^f£SS|L{.= (5.15) 



r2 + {dr/dO) 



In our case, Eq. 1)5. 15(1 



IS 



where 



r Z'"'^'""' {r' - {dr/d9fl (5.17) 
A yjr^ + (dr/d0f 

and r is measured in astronomical units. 

When A > 1/2, because the derivative of the shock with respect to the polar angle 
vanishes for sufficiently large times, the emission flow tends to a constant value given by 



X 



mjyj I 10-7 " i VsOOkms-i/ 



(5.18) 



X — jK , as t ^ oo. 

,104 J V5GHz/ Vl50pc 



Figure lT^ shows the emission flux as a function of time for = 10 ^MQyr ^ (A = 10). 
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§6 Conclusions 

We have developed a simple prescription to describe the hydro dynamical interaction be- 
tween a rotating accretion flow onto a star with a spherically symmetric wind. For simplic- 
ity we considered a purely hydro dynamical flow, with no magnetic fields. The evolution 
of the shock wave produced by the interaction of the wind with the accretion flow is thus 
described by balancing the wind's force ~ M^v^ with the accretion flow's force ~ Mfk, 
where fk is the keplerian velocity that a particle has at a characteristic distance of rj, 
the radius of the disc. That allowed us to deflne a dimensionless parameter A given by 
Eq. 1)3. 5|) . For very large values of A, the surface of interaction expands to infinity and 
the thin layer approximation breaks down. For small values of A, in fact for A < 1/2, it is 
possible to reach a steady configuration. 

We would like to point out that a full analysis of the interaction should take into account 
the gravitational pull that the central star makes on the gas and the equilibrium equation 
must also include a centrifugal term. If magnetic fields are to be included on the whole 
description of the problem as well, then magnetic pressure will also play an important role 
on the structure of the cavities. We will make a more precise analysis considering all these 
physical processes in a future article. 
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Appendix 

Ea. (|2.12j) is a reduced Cardan third order algebraic equation, i.e. 



± 3ax + 2/3 = 0, 



(Al) 



with 



r — 1 



r cos 9 



and X = cos ^o- 



(A2) 



a = 



3 



2 
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I — \r 

The solutions of Eq. ()A2|) are given by f|Namiasl . [l^SJ) 



X = < 



20^1"^ sinh {iarcsinh (— } , 



for a + sign, 



2a^/2 cosh {iarccosh (-/J/a^/s) | ^ 



for a — sign and 0^ — o? > 0, 



(A3) 



2a^/2 cos {i arccos (-/?/a^/^)} , 

for a — sign & — < 0. 
For each case in Eq. HA3|) there exist three different solutions due to the fact that 



arccos z = ± arccosp z + 2A;7r, 
arcsinhz = (— )'^arcsinhpZ + 2A;7ri, 
arccoshz = itarccoshpZ + 2/c7ri, (A4) 

with A; = 0, ±1, ±2, . . . The subindex p refers to the main Riemann sheet. 

For our particular case, we are interested in real solutions, so for the second and third 
cases in Eq. ()A4|) we set A: = 0. For the first case, due to the symmetry of the problem, we 
consider only values of the polar angle such that < < 7r/2. Let us define \ji = —(ijo?!'^ 
and analyse the behaviour of the function C = cos {(arccos /i)/3 + 2A;7r} with A: = 0, ±1. 
When /z > we expect C ^ 0. Since C(At > 0, A; = 0) > then the required solution occurs 
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when k = 0. In other words, the required full solution is given by Eq. with k = 0. 
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Figure 1: Far away from the star, the fluid particles in Uhich's accretion rotate Uke a 
solid body with constant angular velocity (p. These particles follow parabolic trajectories 
near the central object (star). Only upstream trajectories to the equatorial plane are of 
interest. 
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Figure 2: The figure shows the streamhnes of the accretion flow with rotation projected 
on a plane with constant azimuthal angle. Every single streamline is labelled by its initial 
polar angle which is the initial particle's angle with respect to the Z axis as the gas 
begins to be accreted to the star. Lengths in the diagram are measured in units of the 
radius of the disc r^- 
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Figure 3: Accretion density p varies continuously from 9 = io 9 = 7r/2. The density 
tends to infinity in the point r = 1 for the 6 = 7^/2 curve which corresponds to the border 
of the disc. R represents the cyUndrical radius {R = rsin^). Lengths are measured in 
units of the radius of the disc (r^) and the density is measured in units of pQ (see text). 
The polar angle 9' is such that = (7r/60)6''. 
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Figure 4: Density isocontours are shown as continuous lines for the accretion flow. The 
density p is measured in units of the density po- Lengths are measured in units of the 
radius of the disc. Different numbers in the plot represent different values of the density. 



REFERENCES 



35 




Figure 5: The interaction between the stellar wind and the accretion flow produces two 
strong shock waves. It is assumed that both shocks occupy the same position in space. R 
represents the cylindrical radius {R = rsinO), and Z is the polar axis coordinate. Distances 
are measured in units of the radius of the disc ra- 
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Figure 6: The interaction between the stellar wind and the accretion flow forms two shock 
cavities. The figure shows a schematic representation of this interaction with both shocks 
lying close to each other. The arrows in the figure show the direction of the flow velocity 
immediately after crossing the shocks. The cylindrical radius R and the polar coordinate 
Z are both measured in units of the accretion disc radius r^. 
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Figure 7: The post-shock direction is given by the values of the pre-shock tangential 
velocity values for each flow. The tangential velocity of the stellar wind v^t is measured 
in units of the wind's velocity Vy^. The accretion velocity Vr is measured in units of the 
velocity The numbers on the plot are values of the quantity A' defined by A = 0.05A'. 
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Figure 8: The figure shows the shape of the shock for different time intervals t, in units 
of rd/vk and different values of A. The shock converges to the steady case of Figure |21 
when A < 1/2 and for sufficiently large values of time t. When A > 1/2 there is no steady 
configuration and the shock grows indefinitely. The cylindrical radius R and the coordinate 
Z are measured in units of the radius of the disc r^- 
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Figure 9: When A > 1/2, the geometrical form of the shock is spherical for sufficiently 
large times and its radial coordinate grows as a power law (see Eq. p.lll) '). The curves on 
the plot show the shock position r evaluated on the polar axis measured in units of the 
radius of the disc r^. Continuous lines show the numerical calculation and dotted ones 
represent the analytical approximation. The time t is measured in units of rj/vk. 
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Figure 10: When the parameter A > 1/2, the thin layer approximation is no longer valid 
for sufficiently large times. The plot shows the variation of the cooling length (Ic) as a 
function of time. The coordinate r is the distance from the star to the shock. Dotted lines 
correspond to the asymptotic approximation obtained in Eq. (|5.9() and continuous lines 
correspond to the numerical calculation. The time t is measured in years. 
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Figure 11: Emission isocontours normalised to the maximum (shown with crosses). The 
isocontours take values 0.1, 0.2, ... 1.0 with = lO~^M0yr~^ (A = 10) and a resolution 
9b = 0".2 at a distance of 150pc. The angle cl) is formed between the equatorial plane and 
the observer. The time t is measured in years. The cylindrical radius R and the coordinate 
Z are measured in astronomical units. 
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Figure 12: Emission isocontour variations due to different resolutions, normalised to the 
maximum (shown with crosses). The isocontours take values 0.1, 0.2, ... 1.0 with = 
10~^MQyr~^ (A = 10) and t = 3197 yrs. The angle is measured is the angle formed by 
the equatorial plane and the line of sight. The angle 6b represents the resolution. The 
cylindrical radius R and the coordinate Z are measured in astronomical units. 
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Figure 13: When the dimensionless parameter A > 1/2, the emission flux S tends to a 
constant value for t S> 1. The continuous line represents the numerical calculation and the 
dotted hne the asymptotic analytical value. The time t is measured in years and the flux 
S in mJy. 



